Two-gap superconductivity in Lu 2 Fe 3 Si5: a transverse-field muon spin rotation study 



P.K. Biswas, 1 ^ G. Balakrishnan, 1 D.McK. Paul, 1 M.R. Lees, 1 and A.D. Hillier 2 

1 Physics Department, University of Warwick, Coventry, CV4 7AL, United Kingdom 
''ISIS Facility, Rutherford Appleton Laboratory, Chilton, Oxfordshire, 0X11 OQX, U.K. 

(Dated: October 19, 2010) 

The superconducting properties of Lu^FeaSis with T c = 6.1 K have been investigated using 
low-temperature transverse-field muon spin rotation (/iSR) and specific heat measurements. The 
magnetic penetration depth at zero temperature, A(0), is 353(1) nm. However, the temperature 
dependence of the magnetic penetration depth, A (T) is consistent with a two gap s + s-wave model. 
Low-temperature specific heat measurements on the same sample also show evidence of two distinct 
superconducting gaps. 



PACS numbers: 76.75.-K, 74.70.Ad, 74.25. Ha 

The discovery of superconductivity in MgB 2 with a 
T c ~ 39 K— has generated a great deal of interest in su- 
perconducting materials containing light elements such 
as B, C, and Si. Among these materials, the ternary- 
iron silicide superconductors R^FesSis with R = Lu, Y, 
or Sc are particular noteworthy due to the presence of 
ironi^ Lu2FesSi5 is the most interesting of the ternary- 
iron silicide superconductors because of its high super- 
conducting transition temperature (T c = 6.1 K), large 
upper critical field (/j,qH C 2 — 6 T)2^ and unconventional 
superconducting properties. Recently, a detailed study of 
the low-temperature specific heat on a single crystal of 
Lu2FeaSi5 revealed a two-gap superconductivity similar 
to that seen in MgB2.^ 

Lu2Fe3Si5 has a tetragonal Sc2Fe3Sis-type struc- 
ture (space group PA/ rune) consisting of quasi one- 
dimensional iron chains along the c axis and quasi two- 
dimensional iron squares parallel to the basal planed 

Muon spin rotation (/xSR) is a ideal probe with which 
to study the mixed state of type-II superconductors as it 
provides microscopic information of the local field distri- 
bution within the bulk of the sample. It has often been 
used to measure the temperature dependence of the Lon- 
don magnetic penetration depth, A, in the vortex st&teZ& 
The temperature and field dependence of A can provide 
information on the nature of the superconducting gap. 

Here we have investigated the unusual superconducting 
properties of L^FeaSis by carrying out low-temperature 
/iSR measurements on a polycrystalline sample. We show 
that the temperature dependence of A can be well de- 
scribed using a two-gap s + s-wave model. The mag- 
netic penetration depth at T = K is estimated to be 
A (0) = 353(1) nm. We also study the low-temperature 
specific heat of L^FeaSis in order to support the valid- 
ity of the two-gap model. We compare these results with 
published data for the R2FeaSis system. 

A polycrystalline sample of L^FeaSis was prepared 
by melting a stoichiometric mixture of lutetium shot 
(99.99%), iron granules (99.999%) and silicon pieces 
(99.99%) in an arc furnace under an argon atmosphere. 
The as-cast sample was poorly superconducting with a 
T c = 4.8 K and a broad transition. In order to improve 
these characteristics, it is essential to anneal the as-cast 
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FIG. 1: (Color online) The temperature dependence of the 
dc magnetic susceptibility of Lu2Fe3Sis measured using both 
zero-field-cooled warming (ZFCW) and field-cooled cooling 
(FCC). The diamagnetic susceptibility shows a T c onset of 
6.1 K. 



samples at high temperature for a long period of time^^ 
The as-cast sample was sealed in a quartz tube under a 
partial pressure of argon. The sample was then heated at 
a rate of 200° C/h to 800° C, held at this temperature for 
48 h, then heated at the same rate to 1100°C and held at 
this temperature for 72 h. The sample was then cooled 
at 200° C/h to 800° C, maintained at this temperature for 
72 h, and then finally cooled to room temperature, dc 
magnetic susceptibility versus temperature (T) measure- 
ments were performed using a Quantum Design Magnetic 
Property Measurement System (MPMS) magnetometer. 
The temperature dependence of the magnetic suscepti- 
bility shows that the annealed sample of L^FeaSis has a 
transition temperature, T c (onset), of 6.1 K [see Fig. [Q. 
It is the superconducting properties of these annealed 
samples that are discussed below. Low-temperature spe- 
cific heat measurements were carried out using a two-tau 
relaxation method in a Quantum Design Physical Prop- 
erty Measurement System (PPMS) equipped with a 3 He 
insert. 
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FIG. 2: (Color online) The transverse-field muon-time spectra 
(one component) for Lu2Fe3Sis collected (a) at T = 6.5 K and 
(b) at T — 0.3 K in a magnetic field hqH — 30 mT. 



The ^SR experiments were performed on the MuSR 
spectrometer of the ISIS pulsed muon facility. The TF- 
/iSR experiment was conducted with applied fields be- 
tween 5 mT and 60 mT, which ensured the sample was 
in the mixed state. The magnetic field was either ap- 
plied above the superconducting transition and the sam- 
ple then cooled to base temperature (FC), or the sample 
was first cooled to base temperature and then the field 
was applied (ZFC). The MuSR spectrometer comprises 
64 detectors. In software, each detector is normalised for 
the muon decay and reduced to two orthogonal compo- 
nents which are then fitted simultaneously. 

The sample was mounted on a silver plate with a circu- 
lar area of ~ 700 mm 2 and a small amount of diluted GE 
varnish was added to aid thermal contact. The sample 
and mount were then inserted into a Oxford Instruments 
He 3 sorbtion cryostat. Any silver exposed to the muon 
beam gives a non-decaying sine wave. 

TF-/iSR precession signals above and below T c = 6.1 K 
are shown in Figure [21 Above the superconducting tran- 
sition i.e. in the normal state, the signal decays very 
slowly, but the decay is relatively fast in the supercon- 
ducting state due to the inhomogeneous field distribution 
from the flux-line lattice. We can model these inhomo- 
geneous field distributions using an oscillatory decaying 
Gaussian function 



G x (t) = A exp (-At) exp {-a 2 t 2 / 2) cos (wit + <f>) 

+Ai cos (w 2 t + 4>) , (1) 

where U\ and o>2 are the frequencies of the muon preces- 
sion signal and background signal respectively, <j) is the 

initial phase offset, and a is the Gaussian muon spin re- 

i 

laxation rate, a can also be defined as a = (cr 2 c + cr 2 m ) 2 , 
where a sc is the superconducting contribution to the re- 
laxation rate and a nm is the nuclear magnetic dipolar 
contribution which is assumed to be constant over the 
entire temperature range. Fig. [3£i shows the tempera- 
ture dependence of a sc obtained in an applied TF of 0.03 
T. Fig. [3b presents the magnetic field dependence of a 8C 
collected at different temperatures below the supercon- 
ducting transition. A deviation in the field dependence 
of a sc is observed at 40 mT in 0.3 K data. A small devi- 
ation of a sc is also present at the same field in 2 K data, 
whereas it is constant above 2 K. 

The temperature dependence of the London magnetic 
penetration depth, A (T) is coupled with the supercon- 
ducting Gaussian muon-spin depolarization rate, a sc (T) 
by the equation 
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where r y fl /27T = 135.5 MHz/T is the muon gyromagnetic 
ratio and $o = 2.068 x 10" 15 Wb is the magnetic flux 
quantuniili£ A (T) can be calculated within the local Lon- 
don approximatio n 10 i 1 1 by the following expression 
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where / = [1 + exp (E/ksT)] is the Fermi function, 
<p is the angle along the Fermi surface, and A, (T, <p) — 
Ao,j5 [T/T c ) g (ip). The temperature dependence of the 
gap is approximated by the expression 5(T/T C ) — 

tanhjl.82 [1.018 (T c /T - 1)]°' 51 | while g (ip) describes 

the angular dependence of the gap and is replaced by 1 
for both an s-wave and an s + s-wave gap, and |cos (2ip)\ 
for a d-wave gapj 12 ' 13 

The temperature dependence of the penetration depth 
can then be fitted using either a single gap or a two-gap 
model which are structured on the basis of a-mode l 14 ' 15 
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where A~ 2 (0) is the penetration depth at zero- 
temperature, Aq i is the value of the i-th (i = 1 or 2) 
superconducting gap at T — K and is a weighting 
factor with oj\ + ll>2 = 1. 
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FIG. 3: (Color online) (a) The temperature dependence (on 
a log scale) of the superconducting muon spin depolarization 
rate, a sc , collected in an applied magnetic field fioH — 30 mT. 
(b) Superconducting Gaussian depolarization rate, a sc , versus 
applied magnetic field for Lu2Fe3Sis collected below T c at 
0.3 K, 2.0 K, 2.5 K, 3.0 K and 4.0 K. 
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FIG. 4: (Color online) The temperature dependence of the 
London penetration depth as a function of temperature for 
Lu2Fe3Si5. The solid line is a two-gap s + s-wave fit to the 
data while the dashed and dotted lines represent the d-wave 
and s-wave fits respectively. 
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Fits to the data using the three different models are 
shown in Fig. H) The fits appear to rule out the s- 
wave and rf-wave models as possible descriptions for 
Lu 2 Fe3Si5 as the % 2 values for these models are 33.92 
and 15.91 respectively. The two-gap s + s-wave model 
gives a good fit to the data with a \ 2 of 1.94. The 
two-gap s + s-wave model gives Ao^/fceTc = 1.76(6) 
and A . 2 /k B T c = 0.40(1) with u x = 0.35(1). The ra- 
tio of larger to the smaller gap, k, 4.40, which is 
consistent with the value 5 obtained by low-temperature 
specific heat measurement 5 , and 3.44 obtained by pen- 
etration depth measurement using the tunnel-diode res- 
onator technique^ on a single crystal of Lu2Fe3Si5. The 
magnetic penetration depth at T = K is found to be 
A (0) = 353(1) nm. The in-plane penetration depth is 
200 nm, obtained by tunnel-diode resonator technique^ 

We have also performed low-temperature specific heat 
measurements on L^FesSis which support the asser- 
tion that Lu2Fe3Si5 is a two-gap superconductor. Fig. [S] 



FIG. 5: (Color online) The specific heat divided by temper- 
ature (C/T) as a function of T 2 for L^FeaSis. The dashed 
line shows the fit to the data in the normal state. 



shows the specific heat divided by temperature (C/T) as 
a function of T 2 for the same polycrystalline sample of 
Lu2FeaSi5 used for the /xSR study. A pronounced jump 
in the specific heat is observed at 6.1 K which indicates 
that the sample exhibits bulk superconductivity. The 
normal state heat capacity has been fitted up to 12 K 
by C = jT + (3T 3 + aT 5 , where jT is the electronic 
contribution and /3T 3 + aT 5 represents the lattice con- 
tribution to the specific heat. We obtained fitted param- 
eters 7 = 24.9 mJ/mol K 2 , /? = 0.247 mJ/mol K 4 and 
a = 5.38 x 10~ 4 mJ/mol K 6 which are consistent with the 
reported values for both polycrystallineii - — and single 
crystal samples^ We observed a sizeable residual specific 
heat coefficient, 7 = 7.21 mJ/mol K 2 , at T = K. Inter- 
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FIG. 6: (Color online) The temperature dependence of the 
normalized electronic specific heat as a function of T/T c for 
Lu2Fe3Si5. The solid line is a two-gap fit to the data. 



estingly, a finite residual specific heat coefficient has also 
been observed in a polycrystalline sample of the same sys- 
tem^ whereas it is absent in data for a single crystal^. A 
similar effect has also been reported in Bao.6Ko.4Fe2As2 
(70 = 7.7 mJ/mol K 2 ), Ba(Fei_ x C 0a; ) 2 As2 (70 - 3.0 
mJ/mol K 2 ) and Ba(Fe .92Coo.o8)2As2 (70 = 3.7 mJ/mol 
K 2 )i22r— Possible explanations for this residual specific 
heat coefficient involve pair breaking effects of an un- 
conventional superconductor, 23 spin glass behavior, or 
crystallographic defects. 2 - Given the metallurgy of our 
polycrystalline sample and the dramatic effects that an- 
nealing has on the electronic properties, we suggest that 
crystallographic defects are the most likely cause of the 
residual specific heat coefficient in heat capacity data. 

We have calculated the normalized electronic specific 
heat, C e /^T, by subtracting the lattice contribition from 
the total heat capacity and then renormalizing by (7 — 
7o). For more details, see Refs. I22I and l24l Fig. [5] shows 
the temperature dependence of the normalized electronic 
specific heat, C e /jT, for Lu2Fe3Si5 as a function of T/T c . 
We find two clear anomalies in the data which show that 
Lu2FeaSi5 has two energy-gaps. A large jump appears 
at T c and a smaller one at T c /5. The value of C e /jT 
at T c is found to be 1.13(1) meV, which is much smaller 
than the BCS value of 1.43 meV but consistent with the 
value of 1.05 meV measured on a single crystal^ and also 
agrees well with the reported values for polycrystalline 



samples*^— To perform a two-gap fit to the C e /jT data 
in the superconducting state, we use the BCS expressions 
for the normalized entropy, S, and the specific heat 
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j n T c dt 
where t = T/T e , E = [e 2 + A 2 (i)], and y 

The temperature dependence of the energy gap varies 
as A(t) = A S(t), where 5(t) is the normalized BCS 
gapi2£ The solid line in Fig. [5] is a two-gap fit to the 
data. We obtain two distinct superconducting gaps, 
Ai/fc B T c = 2.13(3) and A 2 /k B T c = 0.51(1). The weight- 
ing factor, lj\ — 0.54, which is slightly larger than the 
value obtained from fits to the ^tSR data. The good 
agreement between the experimental heat capacity data 
and the two-gap model argues in favor of the presence 
of two distinct superconducting gaps in L^FeaSis. The 
ratio of the larger to the smaller gap {^) is 4.18, which 
is close to the 4.40 obtained from /iSR measurements on 
the same sample and is also consistent with the published 
data on Lu 2 Fe 3 Si 5 ^^ 

In summary, we have performed a /LtSR study on a 
polycrystalline sample of Lu2Fe3Si5. The temperature 
dependence of the magnetic penetration depth data was 
fitted with three different models. A two-gap s + s-wave 
model provides the best fit to the data. Low-temperature 
specific heat measurements on the same sample also con- 
firm the presence of two distinct superconducting gaps. 
The specific heat results can also be reproduced by a two- 
gap model and support the /xSR results. These results are 
consistent with other reported data for this systemi^i 1 ^— 
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